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/7-mode frequencies in solar-like stars: * 
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Abstract. As a part of an on-going program to explore the signature of p-modes in solar-like stars by means of high-resolution 
absorption line spectroscopy, we have studied four stars (a CMi, n Cas A, £ Her A and /3 Vir). We present here new results from 
two-site observations of Procyon A acquired over twelve nights in 1999. Oscillation frequencies for 1=1 and (or 2) p-modes 
are detected in the power spectra of these Doppler shift measurements. A frequency analysis points out the difficulties of the 
classical asymptotic theory in representing the p-mode spectrum of Procyon A. 
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1. Introduction 

Precise measurements of the frequencies and amplitudes of 
the solar five-minute oscillations have led to detailed infer- 
ences about the solar in ternal structure (cf. review paper by 
IChristensen-Dalsgaarcl2 002 ) . Because of the wealth of knowl- 
edge obtained from helioseismology there is a strong impetus 
to use the pulsation analysis techniques to probe the interiors 
of solar-like stars. From the ground, the application of seismic 
techniques to these stars is difficult because of extremely small 
variations in intensity and velocity associated with p-modes 
and the need to have an adequate temporal coverage to re- 
solve the mo des. To detect su ch small amplitudes, photometric 
methods (Gilliland et alJll993l) suffer from scintill ation noise 
unless conducted from space, e.g. MOST JWalker et alJl2003l 
I Matthews et~ai] l2000h or the future COR OT dBagliiJ l2003l) 
and Eddington ( Roxburgh & Favatal 12003) missions. Recent 
developments in high-precision spectrometric techniques and, 
in particular, Doppler spectroscopy have led to the detec- 
tion of p-mo de fre quencies in ind ividual solar-like sta rs (e.g., 
Martic et a lJll99<* bedding et alJrZOOlf iMartic et alJrZOOlht 



2002; Frandsen < 



lT l2002l iKield: 



sen et al 



Bouchv & Carrier 
2003; Carrier & Bourban 2003). A recen t review of these and 
other measurements has been given by IBedding & Kields'enl 
2003). 
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One of the most interesting target for solar-like seismol- 
ogy is Procyon A (a CMi, HR 2943, HD61421), a F5 IV star 
with ra v =0.363 at a distance of only 3.53 pc. Its fundamen- 
tal stellar parameters are now well determined. Procyon A 
is in a 40-year period visual binary system, sharing the sys- 
temic velocity with a wh ite dwarf. Procyon ' s astrometric orbit 
was recently updated bv lGirard et alJ J2000I) . With a parallax, 
II = 283.2+1.5 mas, they estimated a mass of 1.470±0.045M o , 
which agrees well wit h previous excellent predictio ns from 
stellar interior models (|G uenther & DemaraueJll993l) . In ad- 
dition, iMozurkewich et all dl99 lh . using optical interferome- 
try, measured a stellar angular diameter 9 = 5.5 1+0.5 mas. 
Adopting the very preci sely measured para llax by Hipparcos, 
IT = 285.93 + 0.88 mas. Iprieto et alJ (2002) derived a slightly 
lower mass of 1.42 ± 0.06 M Q , a radius R/R = 2.071 + 0.02 
and a gravity logg = 3.96 ± 0.02 (in cgs units). 

An exploratory sing le-site Doppler observation of Procyon 
A by Ma rtic et alJ ( 119991) . hereafter Paper I, revealed the pres- 
ence of several frequencies that fell in a comb-like pattern with 
roughly equal spacing of 55 pHz in the region of excess power 
around 1 mHz in the power spectrum. Because of the large day- 
time gaps in these single-site measurements, it was not possi- 
ble to identify these frequencies unambiguously. Note however 
that the set of frequencies reported in Paper I were determined 
very close (< 2 u Hz) to the model f requencies published inde- 
pendently by Chab over et al.Nl999l) . 

Following the findings of these exploratory studies we have 
carried out a two-site campaign of Procyon A to identify as 
many p-modes as possible without alias ambiguities. In sec- 
tion |3 we describe the context of the observing runs, in sec- 
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3.00 
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Table 1. Observation log (UT) 



Fig. 1. Doppler shift measurements over seven nights from two- 
sites observations of Procyon A (black dots: ELODIE data, 
green dots: AFOE data). 

tion|3]we present the periodograms of the time sequences from 
coordinated two-site measurements. In sections |4] and [5] we 
outline the different techniques we used to search for the p- 
modes in the power spectra. In the last sections, we present 
the principal mode characteristics (large and small separations, 
amplitudes) and compare the results with the existing models. 

2. Observations and Data reduction 

The observations were carried out with the ELODIE spec- 
trograph mounted on the 1.93m OHP telescope, France and 
the Advanced Fibre-Optic Echelle (AFOE) spectrometer lo- 
cated at the 1.5m (60") Tillinghast telescope of the Whipple 
Observatory, Arizona, USA. The ELODIE and AFOE are fiber- 
fed cross-dispersed echelle spectrographs (Bara nne et al.ll 996: 
iBrown et al.ll 19941) that were optimized for precise and stable 
Doppler measurements. Both systems use double-fiber scram- 
blers. 

As configured for these coordinated measurements, AFOE 
covered 160 nm of spectrum lying between 393 and 665 nm 
(23 echelle orders) at a resolution of roughly R ~ 32000. The 
wavelength reference used for these observations was a ThAr 
emission-line spectrum brought into the spectrograph by a sec- 
ond optica l fiber. The r e duct ion steps for AFOE data are ex- 
plained in Bro wn et al total number of 1127 ex- 
tracted spectra, with sampling rate of about 140 s, has been 
used for this Doppler velocity time sequence analysis. 

The ELODIE data were obtained as a continuous series of 
raw CCD (Tk 1024) frames with a sampling rate of about 100 s. 
We built an asteroseismic mode to automatically acquire long 
uninterrupted sequences of simultaneous stellar and reference 
(Fabry-Perot) exposures. The advantage of using the closely- 
spaced channelled spectrum from a fixed (ZERODUR) Fabry- 
Perot interferometer is that it gives the best possible reference 
even for very short exposure times (30 s for Procyon) and al- 
lows us to monitor with a high precision the spectrograph in- 
stabilities. In our configuration, the reduction process produces 
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Fig. 2. Power spectrum of Procyon A from Doppler shift mea- 
surements over 15 hours, with ELODIE and AFOE. 



for each CCD frame a set of stellar and Fabry-Perot spectra, at 
a resolution of R ~ 42000, interleaved into 67 echelle orders, 
each covering 5.25 nm or in total a wavelength domain from 
390.6 to 681.1 nm. 

The radial velocity compu t ation a lgorithm is based on the 
method proposed by Connes ( 1985) for the computation of 
the very small residual shifts with an absolute accelerome- 
ter. Our algorithm was first tested on i^UMa, HR4335 (see 
IConnes et alJ Fl996) and optimized with a spline interpolation 
for the iterative Doppler shift calculation to take account of the 
larger shifts induced by earth motion and instrument drifts dur- 
ing the night (cf. Paper I). 

The formalism of t he Connes' algorithm was summarized 
■ 1 i 1 ° 

in Bouchv et al. (2001). It was recently u sed as the "optimum 
weight procedure" bv lBouchv & Carried d2.002[) for the detec- 
tion of the oscillations of aCen and bv lFrandsen et alJ (Eo02) 
for the giant star £ Hya. We should however note that the use 
of this method, for close echelle orders, is very sensitive to the 
star flux variation in the course of the night that induces an 
intensity-velocity correlation and in consequence a larger am- 
plitude of the oscillation signal. 
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Fig. 3. Top panels: Power spectra of the ELODIE and AFOE Doppler shift measurements (7 nights) of Procyon. Bottom left 
panel: Power spectrum of the combined data sets. Bottom right panel: Simulated spectrum computed with noise and with the 
same observation times as the combined sequence ELODIE/ AFOE. 



Due to variable weather conditions and lower mean S/N 
than in our 1998 observations, we de-correlate the Doppler ve- 
locity from residual flux variations in and between the orders, 
used for the calculation of the final weighted mean Doppler 
signal. 

A journal of Procyon A observations in 1999 is given in 
Table ^ The table lists the length of the individual runs in 
hours, the number of selected exposures and the rms scatter 
of the time series. To illustrate the precision of each data-set, 
we show in Fig.^the resulting Doppler signal from both sites. 

3. Power spectra analysis 

The power spectrum of the weighted mean Doppler shift mea- 
surements for one of the best nights (high S/N, low rms resid- 
ual) of Procyon observations is presented in Fig.|2] 

Inspection of the other single-night spectra confirms the 
presence of an excess power and several common peaks 
that modulate the envelope between 0.3 mHz and 1.7 mHz. 
Comparison of individual power spectra allows us to search 



for the repetitive occurrence of peaks, free from aliasing, from 
one night to an another, but the frequency resolution is too 
low to distinguish the mode fine structure. We find in a few 
cases large peaks which highly exceed the S/N in the region of 
the hump of power. This characteristic could be explained as 
'large excitat ions' from the studies of the solar p-mode spec- 
trum llElsworth et aljfl995l) . In the solar case, the lifetime of 
the strongest modes are only of the order of a few days and 
their evolution can be followed on a short time scale when their 
frequencies are visibly separated and free from mutual inter- 
ferences. In our observations the occasionally large excitations 
stay visible in the power spectrum from three or four successive 
nights but more or less important interferences with noise have 
also to be considered. In general, the relative amplitudes of the 
frequencies vary from night to night, with the result that a given 
frequency may dominate the amplitude spectrum in some data 
segments and be entirely undetectable in the others. To address 
these issues, we performed many frequency analyses of the en- 
tire data set, numerous subsets of the data and simulated time 
series having the same sampling function as the observations. 
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In Fig. |3 (upper panels), we show the power spectra from 
ELODIE and AFOE one-site measurements of Procyon A over 
seven consecutive nights. Note that the observing run with 
the ELODIE/OHP spectrograph was longer, twelve nights (see 
periodogram in Marti c etallboOlbl) . The power spectra were 
calculated using the Lomb -Scargle (LS) modified algorithm 
jLombl 19761 Scarele 1982) with statistical weights. The mean 
white noise level (cr hf _ ps ), the rms scatter (<x ob ), the frequency 
resolution and computation bin are indicated for each data set. 
The noise level is higher in the AFOE data, which makes less 
evident the excess of power in the low-frequency part of the 
spectra. 

The time series from the two sites were merged by first sub- 
tracting a second order polynomial fit and the clipping data 
over 2 <x for each night. This process yields similar rms scat- 
ter of data from one night to another and allows us to com- 
pute the power spectrum of combined time series without us- 
ing weights (see bottom left panel of Fig. O- To verify that the 
amplitude of the excess signal was not severely reduced by re- 
moving outlying points over 2 <x, we calculated also the power 
spectrum using statistical weights (1/cr 2 ) of the time sequence 
after 4 <x clipping. By comparing these two power spectra, we 
note that the frequencies of the major peaks remain unchanged 
but with more power in the sidelobes of the weighted spectrum. 
For these data the application of the weights has not produced 
an important reduction in noise but rather degraded the window 
function. 

One can see in Fig. [3] that the overall shape of the excess 
power around 1 mHz is not symmetric and does not have the 
typical bell-like envelope of the solar oscillations. This can be 
explained by the effects of gaps/noise and amplitude modu- 
lation that are intrinsic to the star (beating between pairs of 
modes that are closely spaced in frequency and finite lifetime 
of each oscillating mode). The effects of noise interference and 
the window function are visible in the power spectrum (see 
right lower panel of Fig. [5} of the simulated time series with 
an artificial oscillating signal and the rms almost equal to that 
of the observed data sets. 

We performed several simulations in the presence of noise 
with the oscillation signal based on frequencies identified in 
section |5] The amplitude distribut ion of the generated signals 
was produced as in Rarh an et alJ (l999) but with mode (/, m) 
ratios computed for updated Procyon parameters and with the 
maximum mode amplitude of about 35 cm/s near the peak of 
the excess power. 

In the periodogram of the joint ELODIE and AFOE obser- 
vations several isolated peaks are present between 0.3 mHz and 
1 .4 mHz, due to a much better window function (see the ampli- 
tude of the sidelobes in the inset of Fig. 0- These frequencies 
were used as starting ones for the comb power analysis and 
asymptotic fit calculations (cf. next section). 

4. Echelle diagrams 

In the solar case, the asymptotic theory (see e.g.. lTassoulll 980) 
provides quite good results for the acoustic modes with high 
radial order n and low degree I. The asymptotic relation is ex- 
pected to hold for the solar-like p-mode oscillations in other 



stars and gives mode frequencies as 
v(nl) = Av(n + ^/ + ^ + a) + e„i 



(1) 



where Av is the inverse of twice the sound travel time between 
the centre and the surface, a is dependent on the reflective prop- 
erties of the surface and e„i is a small correction term sensitive 
to the central structure of the star. 

In the vicinity of some radial order «o one can parametrize 
the expected mode frequencies as 



v(n p , - v () + Av (n p + ^/) - /(/ + 1)D, 



(2) 



where n p = n - no, vq is the frequency with n = no and / = 0, 
Avo represents the so-called mean large frequency separation 
between p-modes of same degree and adjacent n. Do is re- 
lated to the so-called small frequency separation Sv/j + 2 where 
Svij+2 = v(n, I) - v{n - 1,1 + 2) (4/ + 6)D . 

The third separation cJvoj, which represents the amount by 
which 1=1 modes are offset from the midpoint between the /=0 
modes on either side (see e.g., bedding & Kieldsenll2003l) can 
be defined as <5vo,i = v(n, 0) - j(v(n, 1) + v(n - 1, 1)). If the 
asymptotic relation holds exactly, Do = ISvqz = 4i<? v i3 = 

To the first order, a star pulsating with a number of con- 
secutive radial modes of a given I has a comb-like eigenspec- 
trum with a spacing Avo and if both even and odd I modes are 
present, with spacing Avo/2. To estimate Avq in the region of 
excess power of the Procyon frequency spectrum we used the 
comb response method. 

In previous analysis of the single-site data dMartic et alJ 
2001a) we applied the comb analysis to the CLEAN-ed power 
spectra because of important aliases at spacings ±11 .57 yuHz 
from the largest peaks. To improve the accuracy of these re- 
sults and to see whether we can use the comb response (CR) to 
distinguish the even fro m odd / mode freque ncies, we rewrite 
the function (cf. Eq. 2 in Kieldsen et al. 1995J) as 



C(v ,Av ) = ]~] 



c .Av „„, .Av rN 
S(v + i— +S)S(v -i— +6) 



(3) 



where 6 = ±(z mod 2)Do according to mode degree I — or 
I = 1 at the central frequency vo and a, = 1 for i — 1 or 2, 
a. -, = 0.5 for i = 3 or 4. 

The modified comb responses (Eq.[3} were then calculated 
at the frequencies of the highest peaks in the power spec- 
trum. For each central frequency Vo c , alternatively considered 
as / = and / = 1, we have searched for the maximum CR for 
10 < Av < 80//Hz and 0.3 < D < 2. 

In Fig.|4]we show examples of the modified comb response 
computed for two vo c frequencies at 753 ^Hz and 776 yuHz. 
The responses are scaled to the maximum peak obtained at 
53.5 ^Hz when vo c = 753 fiHz was assumed to be I = 0. 

From the initial set of the highest peaks corresponding to 
the modes, but also to the aliases, the comb analysis gives 
the reduced set of vo c frequencies and corresponding Avo c 
for which the maximum CR is obtained. The upper panel 
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Fig. 6. Example of an echelle diagram from Procyon amplitude spectrum of the joint ELODIE and AFOE radi al velocity mea- 
surements, computed for the modulo mean Avo= 54 //Hz. The frequencies from the standard model (see Cha bover et aljfl999l) 
are indicated by asterisks, plus signs and squares respectively for / = 0,1,2, dashed lines are at / = 0, 1,2+ 11 .6//Hz). The inset 
shows the colour coding applied to the values in each succesive folding frequency range. 
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Fig. 4. Comparison between modified comb responses com- 
puted for two central frequencies of the power spectrum of 
Procyon data (see description in text). 



of Fig. |5] shows the variation of the first-order spacing over 
several different frequency ranges, as determined from comb 
analysis. The arithmetic mean of these spacings is 54 //Hz, 
a value for a large separation in good agree ment w i th the 
theoretical predictions from models JChabover et al.l 1 1999: 
bi Mauro & Christensen-Dalsgaardl200ltlProvost et all2002l) 
In addition, this exploration allowed us, from the comparison 
of the relative power of the best responses, to separate tem- 
porarily the modes by degree. In Fig.|5l(lowe r panel), we show 
these frequencies in the echelle diagram Ccf. lCrec et all 198.3) 
for the modulo mean Avo= 54 //Hz. 

In a second step we calculate the echelle diagram so as to 
preserve the amplitude information. It can be constructed either 
from the power or CLEAN-ed spectra by adding (within cells 
of the order of the frequency bin) the values modulo large split- 
ting, over threshold determined by the mean noise level of a 
time sequence. This means that any points lower than the noise 
level are set at this level. The echelle diagrams constructed in 
this way are used to compare the overall distribution of the 
power for different possible average splittings. In general, the 
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Fig. 5. Top panel: The first-order spacings over several fre- 
quency ranges, as determined from comb analysis. Bottom 
panel: Echelle diagram of the selected frequencies with the best 
comb responses (/ = 0, 1 indicated respectively by diamonds 
and triangles). 



ridge like structure appears for the mean Avo determined from 
previous comb-like analyses. This method was first tested on 
solar data obtained with the same instrument during a day- 
time run. The echelle diagr am of the s i ngle-sit e ELODIE ob- 
se rvations were shown in iMartic et alJ (12001 bl) and discussed 
rn lProvost etall d2002l) . 

We should point out that the distribution of the power in 
these kinds of echelle representations depends greatly on the 
starting folding frequency vo and the threshold level, especially 
when mode amplitudes are decreased by destructive interfer- 
ence in some frequency ranges. In Fig. [6] we show an example 
of such an echelle diagram based on the ELODIE and AFOE 
joint observations. To increase the visibility of the modes in this 
image we modified the power spectrum by a simple normaliza- 
tion to the highest peak in each successive folding frequency 
range. The maximum power is concentrated in the ridge that is 
identified as / = 1 and which shows the offset from 0.5Avo be- 
tween the I - modes. The departure from a linear asymptotic 
relation is evident. It explains why, for the stars like Procyon, 
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Fig. 7. Identification of the /?-modes in the region of the ex- 
cess power of Procyon A. Dashed lines in blue, red and green 
correspond respectively to modes / = 2, 0, 1, 



classical me thods used for mode diagnostics are not easily ap- 
plicable fcf. lBedding et all200ll) . 



5. Mode identification 

The best values, based on our modified comb analysis, for vo c , 
Avq c and Dq c , were then used to reconstruct the parameter- 
ized frequencies in separate frequency regions (vo c ± 2Avo c ). 
These 'model' frequencies have a significant departure from 
the asymptotic relation for one set of mean fit parameters. From 
this step, we have proceeded in an iterative manner to identify 
the modes. 

The observed high amplitude frequencies were classified 
by degree comparing to the model ones and used to obtain new 
fit parameters by performing multidimensional minimization 
of Eq. |2 This method was used because we first selected the 
peaks with the strongest amplitudes, which means that succes- 
sive n p , or I pairs of modes are not always present. We calcu- 
lated the fit for various combinations of mod es with different 
degre e, using the downhill simplex method ( Nelder & Mead 
Il965t) . A certain number of high amplitude frequencies were 
then shifted by ±11 .57 /zHz and included in the selection of the 
modes, which improved the convergence (small max error). In 
Fig. we show the selection of frequencies in the region of 
excess power at 4 cr above the mean white noise, for which we 
obtained the best fit. The values for Avo and Do from the fit 
are : 

Av = 53.6 /uHz 

D = 0.85 //Hz, (5v 02 ~ 6D = 5.1 //Hz 

Given the complexity of this sp ectrum, the noise b ack- 
ground, the short mode lifetimes (cf. Br own et alJ ll991) and 
the possibility of cross-talk between a real frequency and an 
alias etc., we checked the presence of this set of frequencies 
in other independent subsets of the data. We calculated peri- 
odograms from the sequences of three and four consecutive 
nights which have a sufficient resolution (4 //Hz) to resolve 
I — 0, 2 from I = 1 modes. The independent short time se- 
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ries, chosen for their low rms, were processed using CLEAN 
(Rob erts et al.ll987t) and CLEAN-est jFosterl 1995b algorithms 
to minimize the influence of sidelobes in the power spectra. 
The advantage of the second one is that it gives directly an 
initial set of frequencies that can be used for the mode identi- 
fication. Deconvolution methods allow for the identification of 
modes with smaller amplitudes in the spectrum but isolate also 
the aliases instead of the modes. Unfortunately, in the case of 
Procyon the frequency aliases between / = and / = 1 are very 
close to each other, i.e. v/=o + 11.57 =s vi=\ - 11.57. We should 
also note that, whichever subset is used, there are gaps in the 
amplitude spectrum where frequencies in the comb pattern are 
not ex cited to observable am plitudes. 

In iMartic et alJ (1200 1 ah we showed the histogram of the 
major peaks that were considered as recurrent when present 
within ~ 5 //Hz from one sequence to an another using single- 
site data. From a complementary search of the modes in short 
time sequences, we obtained a set of recurrent frequencies for 
further mode analysis of longer sequences where finite life- 
time and random phases from the stochastic nature of mode 
excitation interfere and may reduce the observed amplitudes of 
modes. 

Finally, we explored all the peaks at 2.5 cr above the mean 
white noise level (7cms _1 ) in the amplitude spectrum of the 
time sequence from coordinated two-site measurements. The 
identified modes, with the assigned Z-values (7 = 2, 0, 1 marked 
respectively with squares, asterisks, and plus signs) are pre- 
sented in Fig. [S] with the size of the symbols scaled to the am- 
plitude of the maximum peak in the spectrum. 

For comparison we also show modes (with red symbols) 
from the 1998 observing run (cf. Paper I). Since, the aliases 
are more important in these single site measurements, we used 
the CLEAN algorithm in the data reduction. In some cases, the 
amplitude of the peak is substituted by the higher one of the 
corresponding alias (symbols in circles). The overall distribu- 
tion of the modes confirms and completes the results from two- 
site observations, which increases the confidence in these mode 
frequencies. One can see in Fig.[8]that the / = 1 mode distribu- 
tions are consistent between two sequences while the 1-0,2 
modes are scattered in the range of the 2 //Hz error bar. For a 
given frequency resolution, the uncertainty of the 1-0,2 mode 
identification is due to the decrease of the small separation 6vq2 
with the frequency. 

We list in Table |2 the frequencies of identified modes, de- 
tected more than once in different data sets. Some of the fre- 
quencies, marked in parenthesis, are detected in longer se- 
quences with better resolution but with a more complicated 
window function (see TableQ- In higher frequency ranges, two 
pairs of peaks (1320.8, 1373 and 1328, 1384 //Hz) are recurrent 
in different subsets of the data at the 3 cr level above the noise 
and can be equally identified as / = 1 modes or aliases. Note 
that at the frequencies of the second pair, the modes would fol- 
low an important curvature in the echelle diagram. 

6. />mode analysis 

Different combinations of frequencies of low degree p-modes 
carry information on different parts the of the stellar struc- 
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Fig. 8. Echelle diagram of the identified modes of degrees / = 
0, 1,2 from two-site joint Procyon observations (cf. Table [2). 
For comparison, we indicate with the symbols in red the 
frequencies detected from previous OHP/ELODIE 1998 run. 
Vertical lines correspond to the best fit of the observed fre- 
quencies to the asymptotic relation, with vo=648.4//Hz, Avo = 
53.6 //Hz and D = 0.85 //Hz 

Table 2. Mode frequencies of Procyon A (in //Hz) 
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1 = 1 
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-9 




329. 


352. 
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382. 


(400.8) 403.5 
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427.5 




450.5 (456.6) 


-6 


481.5 


(484.) 


(507.6) 509.7 


-5 


535.5 


541.3 


564.6 (568.) 


-4 


590.8 
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(619.) 623. 


-3 
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648.4 


(671.6) 675. 


-2 
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725.5 (726.3) 


-1 


(749.6) 
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776.7 


nO 
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825.5 or 834.5 


+ 1 


859.1 




884.5 


+2 


(908. 910.) 




941. 


+3 


(964.) 968. 
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+4 


(1021.) 1023.5 


1026.5 


1050.5 (1052.3) 


+5 


(1076.5) 


1080. 


1104.5 


+6 


(1130.) 


1134. 


1158. 


+7 


1186. 


(1190.) 


(1214.5) 


+8 


1241. 


(1244.6) 


1270. 


+9 


1295. 




1320.8 (1328.) 


+ 10 




1353. 


1373. (1384.) 



ture. The large separation between frequencies of modes of 
given degree / and consecutive radial orders is sensitive to the 
global structure and depends on the outer layers of the star. 
It is mainly related to the stellar radius or the dynamical fre- 
quency oc (GM/R 3 ) 5 . Mean separations between frequencies 
of modes with degree, I, I + 2 or I, / + 1, which penetrate dif- 
ferently in the central layers, are mainly sensitive to the struc- 
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Fig. 9. Variation of the large spacing between the modes of 
same degree, compar ed to the predicted cu rve from the stan- 
dard model for I — ( Chabover et al . 1999). The observed av- 
erage large separations for / = 0,1,2 degree are indicated by 
the horizontal lines. 

ture of the deep interior and hence t o the evolutionary state. 
From different models of Procyon by Pro vost et alJ (120021) . in 
the range of considered stellar parameters and overshoot pa- 
rameters, the mean small spacings 5vq2, <5voi, dv^ vary respec- 
tively within 1 //Hz, 1 .5 //Hz and 2 //Hz. Such small differences 
are difficult to measure with the precision of our data. 

In Fig.|9j we show the large separation as a function of fre- 
quency, for I = 0, 1,2. Note that the mean large spacing does 
not depend much on the degree. When the successive orders are 
not detected, the separation was calculated by linear interpola- 
tion. Few mode frequencies are dispersed more than ± 2 //Hz 
around the mean separations. One of the strongest peaks at 
825 ± 2 //Hz, identified as the / — 1 mod e is eventually dis- 
placed by an avoided crossing (see e.g., Audardetal. 1995; 
iDi Mauroetal]l2003l) . but additional data are needed to con- 
firm this possibility. 

Considering the small spacing of about only 4 //Hz above 
0.7 mHz, it is rather difficult with the present data to distinguish 
without ambiguity / = from I = 2 modes of successive radial 
order. We should however note that the number of identified 
I = 2 modes is higher than for I = 0. The observed larger am- 
plitude of these modes is probably due to interferences between 
modes of different m. To estimate the rotational splitting of the 
modes we assumed uniform rotation and the inclination of the 
rotational axis equal to the inclination of the orbital axis of the 
binary system (31.1°) determined by Girar d et alJ ( I2000I) . With 
an estimated radiu s of 2.011 R n and vsiru b etween 2.7km/s 
(Priet o et alJ20 02) and 6.1 km/s Jfiipers the correspond- 
ing rotational splitting is respectively 0.58 and 1 .3 //Hz. We per- 
formed several simulations (cf. section [3} in which synthetic 
oscillation spectra are computed with the rotational frequency 
between these two values. As verified with simulations for a 
given rotational splitting a pair of I — 0, 2 modes may overlap 
and produce one single large peak at shifted frequency in the 
power spectrum. The frequency resolution of joint OHP and 



AFOE data is not sufficient to attempt a determination of rota- 
tional splitting. 

7. Oscillation amplitudes 

Based on Christensen-Dalsgaard & Frandsen (1983) com- 
putations of expected amplitudes from stochastic excitation 
in star s on or near the main sequence, Kieldsen & Bedding 
(1995) proposed a scaling for the oscillation velocity ampli- 
tude y osc oc L/M , where L is the luminosity and M the mass 
of the star. The calculations by Houd eket alJ ( fl9 99) confirmed 
this result, although they fitted better the velocity with V osc oc 
{L/M) s , and s ~ 1.5. However, both the simple scaling law 
and the computation of Houdek et al. grossly overestimate the 
observed amplitudes of Procyon (< 50cm/s). 

To reduce the "exp ected" amplitude of Procyon, 
Kieldsen & Bedding suggested that V osc is inde- 

pendent of the effective temperature amongst stars of a given 
M and radius R. They proposed V/Vq oc go/g (g surface grav- 
ity) which, for Procyon, yields lower velocity ( 70 cm s" 1 ) than 
fiie_scaling V/Vq = 2[(g Q /g) - 6 + (go/g) 4 - 5 ]- 1 JGilliland et alJ 
1993), but still higher than the observations. 

Theoretical amplitudes from models are also overesti- 
mated, ahnost_tlu^e_timeslri^her than the observed one. It ap- 
pears ( Houdek & Gougj] E'o02l) . therefore, that there is some- 
thing wrong with the theory of the pulsations or the convec- 
tion, or their coupling, for this relatively hot and luminous star. 
According to Houdek (2002), for solar-type stars hotter than 
the Sun and with masses M > 1.35 Mq, theoretical damping 
rates may be too small and estimated convective velocities too 
large, leading to predicted amplitudes larger than the data sug- 
gest. 

We note however that the same scaling as in the so- 
lar case, v max ~ 0.6 v ar where t he acoustic cutoff frequency 
Vac 00 g T~g /2 dBrown et alJl99ll) . can be used to obtain a good 
agreement for the frequency at which the oscillation amplitudes 
are largest. Using a recent value for the Trff = 6530 K from 
|PrietoetalJ (l2002). yields v ac ~ 1.7 //Hz and v max ~ 956 //Hz. 

8. Conclusions 

Two-site observations of Procyon A confirmed the presence of 
solar-like oscillations, characterized by a large frequency sep- 
aration Avo = (53.6 + 0.5)//Hz. We identified 45 p-modes 
frequencies in the range 300 - 1400 //Hz with harmonic de- 
grees I = 0,1,2, in good agreem ent with the theoretical pre- 
dictions bv lChabover et alJ )l999tl . The significative departure 
of the modes in the region of the maximum oscillation ampli- 
tudes from the asymptotic relation explains the difficulties in 
finding the mean principal frequency separation Avo and the 
small frequency separation 6vq2 by usual techniques based on 
the best fit of the asymptotic relation from an /-mixed analysis. 
Additional data with higher frequency resolution are needed 
to determine the rotational splitting and the damping time of 
Procyon p-modes. By determining oscillation frequencies with 
an accuracy better than 0.5 //Hz from an extended ground- 
based multi-site campaign on Procyon we will be able to test 
the techniques for the mode detection and prepare the scientific 
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exploitation of future asteroseismic observations from space 
(MOST, COROT, EDDINGTON). 
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